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ABSTRACT 

In this paper we present a photometric and spectroscopic study of the nearby galaxy clus- 
ter Abell 0671 (A671) with 15 intermediate-band filters in the Beijing-Arizona-Taiwan- 
Connecticut (BATC) system and the Sloan Digital Sky Survey (SDSS) data. After a cross- 
identification between the photometric data obtained from the BATC and SDSS, a list of 985 
galaxies down to y ^ 20.0 mag in a view field of 58'x58' is achieved, including 103 spec- 
troscopically confirmed member galaxies. The photometric redshift technique is applied to 
the galaxy sample for further membership determination. After the color-magnitude relation 
is taken into account, 97 galaxies brighter than /ibatc — 19-5 mag are selected as new mem- 
ber galaxies. Based on the enlarged sample of cluster galaxies, spatial distribution, dynamics 
of A671 are investigated. The substructures of A671 are well shown by the sample of bright 
members, but it appears less significant based on the enlarged sample, which is mainly due to 
larger uncertainties in the light-of-sight velocities of the newly-selected faint members. The 
SDSS r-band luminosity function of A671 is flat at faint magnitudes, with the faint end slope 
parameter q:=-1.12. The SDSS spectra allow us to investigate the star formation history of 
bright cluster galaxies, and the galaxies in the core region are found to be older than those 
in the outskirts. No environmental effect is found for metallicities of the early-type galax- 
ies (ETGs). Both mean stellar ages and metallicities in bright member galaxies are found to 
be correlated strongly with their stellar masses assembled, and such correlations are depen- 
dent upon morphology. The possitive correlation between age and stellar mass supports the 
downsizing scenario. By comparing ETG absorption-line indices with the state-of-art stellar 
population models, we derive the relevant parameters of simple stellar population (such as 
age, [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe]). The ETGs at cluster center tend to have 
smaller H/3 indices, indicating that central ETGs are likely to be older The distribution of 
total metallicity indicator, [MgFe]', does not show any environmental effects. The relations 
between the simple stellar population parameters and velocity dispersion in A67 1 are in good 
agreement with previous studies. 

Key words: galaxies: clusters: individual (A67 1) - galaxies: distance and redshifts - galaxies: 
kinematics and dynamics - galaxies: evolution - methods: data analysis 



1 INTRODUCTION 

According to the hierarchical scenario of structure formation, mas- 
sive clusters form by merging small grou ps continuously and ac- 
creting field galaxie s along the filament ( I West et al]|l99ll , Il995l : 
IColberg et alTlaOOOl) . Optical cluster surveys reveal that many 

i ;alaxy clusters have evidence for dy namically bound substructures 
Rhee et al.ll99ll : lBeers et aljl99lll . A significant fraction (~40%- 
50%) of clusters show multiple peaks or irregular surface bright- 
ness distribution in the X-ray images, indicating that they are still 
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at dynamically active stage, far from equilibrium jjones & FormanI 
[1999; Schuec ker et al . 2001). Compared with the Einstein-de Sit- 
ter case, clusters in early-epoch universe are expected to be more 
relaxed and le ss substructured, as supported by many N-body sim- 
ulation works jcrone et al]|l996l ; lThomas et al.lll998t) . The fraction 
of substructured clusters at different redshifts is thus a useful statis- 
tical quantity directly relevant to cosmology. Studies on the dynam- 
ics of galaxy clusters thus provide a unique tool to put constraint on 
the models of cluster formation and evolution. 

Dense environment in galaxy clusters should have produced 
influence on physical properties and evolutionary path for the mem- 
ber galaxies. Previous studies have found that the observational 
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propertie s of galaxies correlate strongly with local galaxy envi- 
ronment dOomez et al1l2003l : iKauffmann eTZIl2004l : Isaldrv et alj 
l200d h. One of the most well-studied relation in galaxy cluste rs is 
the morphology-density relation (Dressier 1980 ; Post man & Gellet 



1984 IWhitmore & Gilmore|[T99TI ; looto et aljboosi : iHolden et al 



20071) . The core region of a cluster is usually dominated by early- 
type galaxies (ETGs), while the outer region is dominated by late- 
type galaxies (LTGs). It is well appreciated that the LTGs gradu- 
ally lost their gas reservoirs when they were accreted into the core 
region, and finally evolved into lenticular galaxies (SO). This pic- 
ture of morphology evolution in galaxy clusters has been supported 
in high-z morphology-density relation studies fcr essler et al. 1997; 
iFasano et al.ll2000l ; ISmith et alj|2005l ; IPostman et al.ll2005 ). How- 
ever, it is still uncertain and controversial how local galaxy envi- 
ronment affects star formation histories of cluster galaxies. 

The Beijing-Arizona-Taiwan-Connecticut (BATC) system has 
spent much time in observing a sample of more than 30 nearby 
(z < 0.1) galaxy clusters at different dynamic statuses, aiming 
at studying their dynamic substructures, luminosity functions, and 
the star formation properties of cluster galaxies. Abell 0671 (A671; 
z=0.0502) is one target of the BATC galaxy cluster survey. Its Abell 
richness R is set to be ( Abell 1958), with Bautz-Mogan type II-III 
( lBautz&Morganlll970l) . The X-ray emission from the cluster cen- 
ter has been detected by the Einstein observatory and the ROSAT 
All-Sky Survey (RASS). The X-ray luminosity of A671 detected 
in the RASS 0.1-2.4 keV ba rid is 0.9 x l O'^'^erg s"\ and the X- 
ray temperature is 3.1 keV dEbeling et al.lll998 ). which confirms 
that this cluster is a relatively poor system. Figure 1 shows the 
smoothed contours of the Einstein X-ray image and the radio map 
at 1.4 GHz from the NRAO VLA Sky Survey (NVSS), superim- 
posed on the optical image in the BATC— /i band. No radio emis- 
sion is detected at the center of A671. The X-ray surface density 
co ntour is quite regular w ith a single synmietric peak as reported 
by I Jones & FormM jl999l) . However, the detail structure of A671 
is possibly blurred due to the low resolution (~ 1 arcmin spatial 
resolution) and the large PSF (FWHM ~ 1 .5 armin) of Einstein IPC 
image . It is easily seen that the X-ray emission peak does not co- 
incide with the central brightest galaxy, IC 2378, with a positional 
offset of about 90 kpc. A67I is included in the cluster sample of 
the Wide-field Nearby Galaxy-cluster Su rvey (WINGS ) , and two 
substructures in A67I have been found bv lRamella et al.l ( 20071) re- 
cently. For a better understanding of the dynamics of A67I, it is 
important to construct a large sample of member galaxies, and the 
faint galaxies (18.0 < ttIv < 19.5) should be taken into account. 
In this paper, we present a multicolor photometry of the galaxies 
in A67I region with the Beijing- Arizona-Taiwan-Connecticut sys- 
tem. We try to enlarge the sample of cluster galaxies by applying 
the photometric redshift technique to the spectral energy distribu- 
tions (SEDs) of the BATC-detected faint galaxies. Based on the 
SDSS spectra of bright member galaxies, we will derive the star 
formation histories and chemical abundances, and try to find any 
clues of environmental effects on the physical parameters of clus- 
ter galaxies. 

This paper is organized as follows: we present the BATC pho- 
tometric observations and data reduction in Section 2. In Section 3, 
we analyze the galaxies with known spectroscopic redshifts in the 
A67I field. In Section 4, we apply photometric redshift technique 
to select faint member galaxies in A67 1 . In Section 5, dynamic sub- 
structures and luminosity function are investigated based on the 
enlarged sample of member galaxies. In Section 6, we derive the 
star formation histories and chemical abundances of the ETGs in 
A67I. Finally, we summarize our work in Section 7. Throughout 
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Figure 1. The smoothed contours of Einstein image (0.5 — 4.5 keV 
band) (red line) and the NVSS map at 1.4 GHz (green line), superimposed 
on the BATC— h band image. The sizes of gaussian smoothing windows ai'e 
adopted as 30" and 1.2' for radio and X-ray contours, respectively. 



this paper, we assume the cosmological parameters as Qn-i ~ 0.3, 
= 0.7, ffo = 70 km s-^Mpc'V 



2 OBSERVATION AND DATA REDUCTION 

The BATC survey is based on the 60/90 cm //3 Schmidt Tele- 
scope of National Astronomical Observatories, Chinese Academy 
of Science (NAOC), located at Xinglong Station. The BATC sys- 
tem contains 15 intermediate-band filters, covering a wavelength 
from 3000 to 10000 A, which are designed t o avoid night sky 
emission lines dFan et alJl993;lKong et alj200(]|). The transmission 



curves of BATC filters can be found in iFan et alj ( Il996l) . Before 
October 2006, a Ford CCD camera with a format of 2048 x 2048 
was mounted on the telescope, and photometric observations in 12 
bands, from d to p, were carried out. The viewing field was about 
58' X 58', with a scale of l".7/pixel. For pursuing better spatial 
resolution and higher sensitivity in three blue bands, a — c, a new 
E2V CCD with 4096 x 4096 pixels was then equipped. The field of 
view becomes larger (92' x 92') with a spatial scale of l".35/pixel. 
The newly equipped CCD camera has a high quantum efficiency of 
92.2%. 

From March 2003 to October 2007, we totally accumulated 50 
hours exposure for A671 with 15 filters (see the observational in- 
formation in Table 1). With a n automatic data-processing software, 
PIPELINE I ( lFanetalJl99^ . we carry out the standard procedures 
of bias subtraction, flat-field correction, and position calibration. 
The technique of integral pixel shifting has been used in the image 
combination during which cosmic rays and bad pixels are removed 
by comparing multiple images. 

For detecting and measuring the flux of the sources within a 
given aperture in the BATC images, we convert the a — c com- 
bined images to make the pixel size identical with d — p im- 
ages. We adopt a radius of 4 pixels as a photometric aperture for 
all the BATC images to the s ources detected by the SExtractor 
codes dBertin & Amoutslll996h . The flux calibration of SE Ds is 
performed by using the Oke-Gunn dCunn & Strvkeilll983l) Stan- 
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Table 1. The detail of the BATC filters and observation information of A67 1 



Number 


Filter 
name 


(A) 


FWHM 

(A) 


Exposure 
(second) 


Number of 
Images 


Seeing" 
(arcsec) 


Objects 
Detected 


Limiting 
mag 


1 


a 


3369 


222 


18000 


15 


3.79 


3945 


21.5 


2 


b 


3921 


291 


7200 


6 


4.90 


4158 


21.0 




c 


4zUj 


3oy 


12600 




A /in 


5738 


O 1 o 

21. U 


4 


d 


4550 


332 


19500 


17 


4.61 


4739 


20.5 


5 


e 


4920 


374 


15000 


14 


3.78 


5453 


20.5 


6 


f 


5270 


344 


13800 


13 


4.40 


5466 


20.0 


7 


g 


5795 


289 


7800 


8 


4.12 


5506 


20.0 


8 


h 


6075 


308 


7500 


7 


5.62 


5723 


20.0 


9 


i 


6660 


491 


5460 


7 


3.51 


6675 


20.0 


10 


j 


7050 


238 


7500 


7 


4.14 


6017 


19.5 


11 


k 


7490 


192 


12600 


12 


4.09 


5968 


19.5 


12 


m 


8020 


255 


11100 


10 


4.74 


6176 


19.0 


13 


n 


8480 


167 


11100 


10 


4.24 


5893 


19.0 


14 





9190 


247 


15000 


14 


4.13 


5846 


18.5 


15 


P 


9745 


275 


18600 


16 


4.46 


5046 


18.5 



" seeing of the combined image 



dard stars (HD19445, HD84937, BD+26d2606 anA BD+17d4708) 
which were observed during photometric nights. The detaile d in- 
formation about cahbration can be found in lZhou et alj j200lh . Be- 
cause we have no calibration images for a, b, c filters, we instead 
perform the model calibration that has been dev eloped specially for 
the large-field photometric system bv IZhou eTa l. ( 1999). As a re- 
sult, the SEDs of 6782 sources have been obtained in our catalog. 
By cross-identifying the BATC sources with the SDSS photomet- 
ric data within a search circle (defined as a circle of 1.5 arcsec), all 
sources are classified into galaxies and stars. As a result, 985 galax- 
ies brighter than the BATC— /i band magnitude limit are found by 
both surveys, which offers a sample for further analysis. 



3 ANALYSIS OF GALAXIES WITH KNOWN 
SPECTROSCOPIC REDSHIFTS 

3.1 Distribution of spectroscopic redshifts 

For studying the dynamics of galaxy cluster A671, 205 galaxies 
with known redshifts in our viewing field are extracted from the 
SDSS DR8 galaxy catalog. Figure 2 shows the distribution of spec- 
troscopic redshifts of these galaxies. The main concentration with 
a peak at 2 ~ 0.05 is isolated and less contaminated. There are 103 
galaxies with 0.04 < z < 0.06 and they are selected as the member 
galaxies of A671, to which we refer as Sample I. To characterize the 
velocity distribution, we convert the spectroscopic redshifts (zsp) 
into the rest-frame velocities (u) by u = c x (zsp — Zc)/(1 + Zc), 
where c is the light speed, Zc is the cluster redshift with respect to 
the cosmic background radiation. We take the NED-given cluster 
redshift Zc =0.0502 for A671. The velocity distribution can be fit- 
ted by a Gaussian with a dispersion of a = 625 km s^^, and it is 
significantly deviated from a standard Gaussian function. 

To qualify the distribution of rad ial velocities of m ember 
galaxies, we use the ROSTAT software l lBeers et al.|[l990h to cal- 
culate two resistant and robust estimators, biweight location (Cbi) 
and scale (Sbi), analogous to the velocity mean and the standard 
deviation. For these 103 galaxies, we achieve Cbi = 14561^32" 



km s-\ Sbi = 820^50 km s'V The errors are determined from 
the 68% confidence limits based on 10,000 bootstrp resamplings 
of the velocity data. Taking a cosmological correction factor of 
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Figure 2. Distribution of the redshifts for 205 galaxies in the A67 1 field, 
with bin size Az = 0.001. The galaxies in A671 are centered at z ~ 
0.05. We select the galaxies with 0.04 < z < 0.06 as member galaxies. 
The smaller panel shows the distribution of rest-frame velocities of member 
galaxies in detail, and dash line represents a gaussian fit to the histogram. 



{1 + z) into acc ount, the velo c ity di spersion of A671 should 
be 7801:^7 km s~ MAguerri et al.1 i2007h studied a sample of 88 
nearby clusters. Only 72 galaxies were included in their sample, 
and a smaller Sbi was derived. They found A671 having Cbi = 
14599^33 kms"^and Sbi = 6Wt.ll km s"\ Our statistics is 
based on a larger sample, thus is more reliable. Both studies con- 
firm that A671 has a comparatively small velocity dispersion. 
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Figure 3. Left: Spatial distribution of 103 known member galaxies of A671, including 63 eai'ly-type galaxies (red circles), and 40 late-type galaxies (blue 
triangles). The contour map of surface density for all these galaxies, smoothed by a Gaussian window with a = 1.6 arcmin. The contour levels are 
0.09,0.15,0.21,0.27,0.33 and 0.39 arcmin^^, respectively. Right: Bubble plot for groups of six nearest neighbors, showing the localized variation in velocity 
distribution. 
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Figure 4. Stripe density plot of velocities of the spectroscopically con- 
firmed galaxies in whole cluster, clumps A, B, and C, respectively. 

3.2 Spatial distribution and localized velocity structure 

Because A671 is a nearby cluster, our BATC viewing field can not 
cover whole cluster region. Our photometry focus on a central field 
of 3.4 X 3.4 Mpc'^. The radius of galaxy cluster, r2oo, is defined 
in former studies as the boundary of a cluster, within which the 
mean inn er density i s 200pc, where pc is the critical density of the 
Universe l lGotlll972h . We calculate the r2oo for A671 following the 
formula sues ested bv lCarlberg etai] l ll997l) . The r2oo is a function 
of velocity dispersion. By applying the Sm we derived, the r2oo 
of A671 is 2.25 Mpc, corresponding a slightly larger area than our 
viewing field. 

Before studying the dynamic structure of A671, we try to clas- 
sify the 103 known members into ETGs and LTGs. The early- 
type members should meet two requirements: (l)with no evident 
emission lines;(2)with no evident galaxy arms. We firstly extract 



those galaxies with EW(Ha)< 5 A as early-type candidates. The 
EW(Hq) values are taken from the MPA/JHU catalo£] of SDSS 
galaxies. Then we inspect their SDSS-given images and removing 
those with arms. As a result, the 103 galaxies are classified into 63 
ETGs and 40 LTGs. 

The left panel of Figure 3 presents the spatial distribution of 
103 known member galaxies within our field of view, with the cen- 
tral position of A671, R.A.=8''28'"29', Dec.=30°25'0l" (for the 
J2000 equinox). We superpose the contour map of surface density 
that has been smoothed by a Gaussian window with a = 1.6'. As 
shown in Figure 3, the member galaxies are mainly concentrated in 
the central region within a radius of 1 Mpc. The irregular contour 
in the east and north corresponds to the two substructures found by 
iRamella et al] ( |2007|) . More statistical tests should be performed 
before we can reach a firm conclusion that there are substructures 
in A671. 

To show the substructures of A67 1 in both velocity space and 
projected map, we make use of the K-test CColless & Dunn 199^ 
for the 103 galaxies. The statistic variable, Kn, is defined to quantify 
the local deviation on the scale of groups of n nearest neighbors. A 
larger Kn indicates a greater probability that the local velocity dis- 
tribution differs from overall velocity distribution. The probability 
P{Kn > can be calculated by Monte Carlo simulations with 

random shuffling velocities. When the scale of the nearest neigh- 
bors n varies from 3 to 9, the probabilities P{Kn > k"^^) are 
nearly zero, which means the substructure appears very obvious 
at different scales. The bubble plot at the scale of n = 6 is given 
in the right panel of Figure 3. Since the bubble size is proportional 
to — logP(Z)n > Df-^^), the clustering of large bubbles is a good 
tracer of dynamical substructure. 

As we can see in the bubble plot, the central region of A671 is 
dominated by two clumps of bubbles, and a clump of large bubbles 

^ http://www.mpa-garching.mpg.de/SDSS/DR7 
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in the north-east is also remarkable. We refer to these three clumps 
as A, B, and C, which contain 16, 10, 8 galaxies, respectively. For 
confirming whether these clumps trace the real substructures, we 
present the velocity distributions for subsamples A, B, and C in 
Figure 4, as well as the velocity distribution of the whole sample. 
It is easily seen that the velocity distributions for the three clumps 
are indeed deviated from that of the whole sample. The mean ve- 
locities of the subsamples A, B, and C are 14638lg4i, 15567i26e, 
and 13992lJ^4 km s^^, respectively. The mean velocities of sub- 
samples B and C are significantly deviated from the mean veloc- 
ity of the whole sample, 14561 km s~^. Though the mean velocity 
in clump A is similar to that of the whole sample, the galaxies in 
clump A have a remarkable bimodal velocity distribution, which 
may imply the existence of two different groups. To quantify the 
bimodality signi ficance, we apply the Gau ssian mixture modeling 
(GMM) method jMuratov & Gnedinll2010h to the velocity distribu- 
tion for clump A. The result shows that clump A consists of two 
components, which peaked at 13829 km s^^and 15824 kms^^. 
The unimodal velocity distribution can be rejected at 99% signifi- 
cance level. 

Although the velocity distribution of the subsamples are sig- 
nificantly deviated from the whole sample, their masses should also 
be large enough if they are real substructures. The masses of the 
A671 and its clumps can be estimated by applying the virial theo- 
rem. Assuming that each subcluster is bound and the galaxy orbits 
are random, the virial mass (A/t,) can be derived from the following 
standard formula (Gelle r & Peebles 1973 ; .Oegerle & Hilll994.) : 

N 

where Ur is the light-of-sight velocity dispersion, D is the cosmo- 
logical distance of the cluster, Np — N{N — l)/2 is the number 
of galaxy pairs, and 9ij is the angular distance between galaxies 
i and j. The viral masses of A671 is 9.49 x lO^M©. Clump B 
is the most remarkable among the 3 subsamples, with viral mass 
of 1.53 X The A clump consists of two groups, with vi- 

ral mass of ~ 7.0 x lO^^M© for each. Clump C has similar viral 
mass. Thus we conclude that A671 is not a simple relaxed cluster, 
but most likely at a dynamically active stage. 



4 SED SELECTION OF FAINT CLUSTER GALAXIES 

The technique of photometric redshift can be used to estimate 
the redshifts of galaxies by using the SED information covering 
a wide range of wavelength instead of spectroscopy. This tech- 
nique has been extensively applied to the multicolor photometric 
surveys for detecting the faint and dis tant galaxies, and for subse- 
quential selection of clus ter galaxies jFemandez-Soto et alj|l99^ : 
lllbert et alj|2009l : iKong e t al. 2009) . Based on the standard SED - 
fitting code called Hyper Z (Bolz onella. Miralles & Pelloll2000h . 
for a given object, the photometric redshift, Zph, corresponds to the 
best fit (in the -sense) of its photometric SED with the template 
SED generated by convolving the galaxy spectra in the template 
library with the transmission curves of specified filters. Previous 
work has demonstrated the accuracy of photometric reds hift with 
the BATC multicolor photometric data ( X ia et al. 2002; Ya ng et al] 
|2004 iLiu et ai]|201 ihlzhang et al.ll201 11) . In our SED fitting, only 
normal galaxies are taken into account in the reference tem plates. 
Dust extinction with a reddening law of the Milky Way jAUenl 
Il976h is adopted, and Av is allowed to be flexible in a range from 



0.0 to 0.5, with a step of 0.05. The photometric redshift of a given 
galaxy is searched from 0.0 to 0.6, with a step of 0.005. We apply 
this technique to all BATC galaxies brighter than ftsATC = 19"* .5. 
The procedure of SED-fitting has provided the best-fit photometric 
redshift and its uncertainty for each galaxy. 

From the 205 galaxies with known spectroscopic redshifts 
(zsp), we select 167 galaxies (including 91 member galaxies) that 
are simultaneously detected in at least 12 BATC bands to derive 
their values. A comparison between and z^p values is 
shown in the top-left panel of Figure 5. The error bar of Zph corre- 
sponds to 68% confidence level in photometric redshift determina- 
tion. The solid line denotes Zph = Zsp, and the dashed lines show 
an average redshift deviation of 0.02(1 + z). It is obvious that our 
Zph estimates are basically consistent with their z^p values. For 91 
member galaxies in our spectroscopic sample, the mean value and 
standard deviation of their Zph values are 0.0505 and 0.0112, re- 
spectively. There is no systematic offset in the Zph domain with 
respective to the z^p distribution. Statistically, 85 member galax- 
ies (about 93 percent) are found to have their photometric redshifts 
within ±2(j deviation, in a range from 0.028 to 0.073, demonstrat- 
ing the robustness of our Zph estimate. This Zph region can be ap- 
plied as a selection criterion in the following membership determi- 
nation for faint galaxies. 

The top-right panel of Figure 5 shows the histogram of pho- 
tometric redshifts for all galaxies down to /iBATC = 19™. 5. The 
black histogram shows Zph distribution for the 167 galaxies with 
known spectroscopic redshifts mentioned above. As expected, the 
peak in the Zph distribution is around z=0.05. In the bottom-left 
panel of Figure 5, we show the Zph uncertainties for these 167 
galaxies as a function of BATC— band magnitude. It is remark- 
able that the Zph deviation of fainter galaxies tends to be larger. 
For the faint galaxies with /ibatc ~ 18™. 0, our Zph estimate is 
still robust, but with larger uncertainty. For the remainder galax- 
ies without Zsp values, we give a plot of their Zph uncertainties 
versus BATC— band magnitudes in the bottom-right panel in Fig- 
ure 5. The larger Zph uncertainties for faint galaxies are mainly 
due to larger magnitude errors in photometry. For a reliable mem- 
bership determination based on the Zph estimate, we exclude the 
galaxies fainter than /ibatc ~ 19™. 5, and take the galaxies with 
0.028 < 2ph < 0.073 as member candidates. Due to the robustness 
of our photometric redshift technique, it is conservative that our se- 
lecting criterion would be able to select 80-90% of faint members 
with least contaminants. 

It is well known that there exits a correlation betwe en color 
and a bsolute magnitude for the ETGs (C-M relation) ( Bower et al.l 
Il992l) . in the sense that brighter ETGs appear redder, which can 
be used for verifying the membership selection of the ETGs. The 
left panel of Figure 6 presents the correlation between the color 
index, b ~ h, and BATC— /i band magnitude for the whole mem- 
ber candidates down to /ibatc = 19™. 5, while the right panel 
shows the SDSS color index, g — r, and r— band magnitude. The 
diagrams include the following categories of sources : (1) spec- 
troscopically confirmed early-type member candidates (denoted by 
filled circles), (2) spectroscopically confirmed late-type member 
candidates (denoted by filled triangles), (3) newly-selected faint 
member candidates (denoted by open triangles). The solid line rep- 
resents the linear fitting with the 63 spectroscopically confirmed 
ETGs: b~ h = -0.09(±0.02)/i + 3.47(±0.43), and the dashed 
line represents la deviation. The linear fit of the ETGs in the right 
panel of Figure 6 is g - r = -0.026(±0.006)r- + 1.21(±0.10). 
As shown in both figures, the early-type member galaxies follow a 
very tight color-magnitude relation, faint member candidates also 
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Figure 5. Top-left: comparison between Zph and Zgp 167 galaxies detected with at least 12 BATC bands. The error bars shows the 68% significance coefficient 
level. Solid line denotes Zph = %p, and the dash lines denote the deviation of 0.02(1 + z). There is no systematic deviation between the photometric redshifts 
and the spectroscopic redshifts. Top-right: Zph distribution for the galaxies down to ^batc = 19"' -5, with a bin size of = 0.004. The black histogram 
shows the Zph of the 167 galaxies with Zsp. The dash line shows our faint member candidates selection criterion. Bottom-left: difference of Zph and Zap for 
the 167 galaxies with spectroscopic redshifts. Circles denote member-galaxies, and triangles denote non-members. The dash lines denote the member selection 
cataria(2<T = 0.023). The error bars are given by the Hyper Z codes, with 68% confidence level. Bottom-right: crosses denote la errors of the of the Zp\i for 
the left galaxies with h-band magnitude brighter than 19.5 mag in the viewfield of A671. 



foUow the same C-M relation basically, but seem to be more scat- 
tered. This might be caused by some high-z galaxies that have been 
mixed into our faint member candidates. For excluding these con- 
taminants, we utilize the SDSS C-M relation, and remove the ETG 
candidates with color indices g — r 0.15 mag redder than the red 
sequence denoted by the black solid line. 

Finally, we obtain a list of 97 newly-selected member galax- 
ies. Combined with those 103 spectroscopically confirmed mem- 
bers, we form an enlarged sample of 200 galaxies in A671, to which 
we refer as Sample II in following investigation. 



5 THE PHYSICAL PROPERTIES OF A671 

5.1 The Spatial Distribution and Velocity Structure 

The left panel of Figure 7 shows the projected positions of the 
galaxies in Sample n, superposed with the contour map of sur- 
face density smoothed by a Gaussian window with a = 1.6'. The 



103 member galaxies with known spectroscopic redshifts are de- 
noted by filled symbols, and 97 photometrically selected galaxies 
are denoted by open symbols. Red symbols represent the early- 
type member galaxies and blue symbols represent late-type ones. 
Our BATC multicolor photometry facilitates the finding of large 
number of faint member galaxies, and makes the underlying sub- 
structure along north-east direction more remarkable. Basically, the 
distribution of faint galaxies traces that of bright ones, and no sig- 
nificant substructures are found with Sample II. 

The morphology segregation becomes more remarkable in 
Sample II. Both bright and faint ETGs are highly concentrated 
in the core region, while the LTGs are scattered in the outskirts. 
The shape of contour map seems in accord with the X-ray image, 
which might demonstrate the reliability of our membership selec- 
tion based on the BATC multicolor photometry. 

For detecting the potential substructures in A67 1 , we perform 
the K-test for Sample II. The right panel of Figure 7 shows the 
bubble plot that characterizes the degree of difference between the 
localized velocity distribution, for groups of six nearest neighbors. 
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Figure 6. The left panel shows the BATC C-M relation, using color index b — h and /i— band magnitude, for the spectroscopically confirmed members and 
the galaxies with 0.028 < 2ph < 0.073. The right panel shows the SDSS C-M relation, using color index g — r and SDSS r— band model magnitude. 
The filled red circles denote the spectroscopically confirmed ETGs in A67 1 . The filled triangles denote the spectroscopically confirmed LTGs in A67 1 . The 
open triangles denote the newly selected galaxies with 0.028 < Zph < 0.073. The solid red lines represent the subsequence fitted with the known early-type 
members. The dashed red lines represent the Icr of the best-fit C-M relations. The black solid line represents our selection criterium based on the SDSS C-M 
relation. 
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Figure 7. Left: spatial distribution of all member galaxies in Sample 11. Filled symbols denote the galaxies with known spectroscopic redshifts, while open 
symbols represent newly selected members. Red ckcles denote red-sequence galaxies, while blue open triangles denote those bellow Icr of the red-sequence. 
Right: bubble map shows the local velocity distribution for groups of six nearest neighbors for the whole member galaxies. 



and the overall velocity distribution. We performed 10^ simulations 
to estimate probability P(«:n > «:°'"^) for different group sizes. The 
probability is found to be more than 5% in all cases, which means 
that no substructure is detected at 2a significance (see Table 2). 
This is not consistent with the conclusion that we have achieved 
based on Sample I. We think that the substructure unveiled by the 
spectroscopic redshifts is surely true. Above inconsistency can be 
well interpreted that the velocities derived from the Zpi^ estimates 
are not accurate enough to reflect the subtle velocity structure. The 
abnormity in velocity distribution of substructures might have been 



smoothed/swept by the Zph uncertainties of 97 newly selected faint 
galaxies. So the K-test on Sample II might be misleading. Follow- 
up spectroscopy of these faint member galaxies are needed if one 
wants to investigate the substructures in A67I in detail. 



5.2 The Luminosity Function 

The luminosity function (LF) is a key diagnostic for clusters be- 
cause it is tightly related to dynamical evolution and merging his- 
tory of galaxy clusters. The LF has been widely studied in the 
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Table 2. Result of the K-test for member galaxies 
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Figure 8. Top left: Density map of a 12 X 12 Mpc region around A671. 
Contour levels ai'e from Ict to 5cr, spaced Icr apart. Top right: The region 
which is used to estimate the local background. Dots denote the galaxies 
with the SDSS r-band magnitude brighter than 21.0 mag. Empty regions 
represent the removed overdensity area. Bottom left: galaxy number counts 
extracted within r=1.7 Mpc region (filled circles), and those in the local 
background (open circles) for A671. Bottom right: background-subtracted 
number counts for A671. The solid line represents the best-fit Schechter 
function. 



lect a set of cluster member candidates by utilizing the photomet- 
ric redshift technique. As seen in Figure 5, the accuracy of Zph is 
a function of galaxy apparent magnitude. Further corrections are 
still needed to remove the contribution of contaminant sources and 
compensate the missing members when investigating the LFs. Un- 
fortunately the exact form of the correction function, particularly at 
the faint magnitudes, is very difficult to be derived. 

Using the SDSS r-band photometric data, we perform the sta- 
tistical background subtraction to estimate the contribution of non- 
members to the number counts of galaxies in the cluster direction, 
by measuring the projected number counts of field galaxies outside 
the cluster region. The background is estimated with a 12 x 12 
Mpc region centered on the cluster centroid , outside the clus- 
ter region defined by a radius of 3 Mpc, where the contamination 
fr om cluster galaxies s hould be negligible. Following the method 
of lPaolillo et al.l ( 1200 ll) . we firstly generate a density map of galax- 
ies in the background region by convolving the projected distribu- 
tion of galaxies with a Gaussian kernel of (j=250 kpc in the cluster 
rest frame (the typical size of a cluster core). Then we mask out 
all density peaks which are above 3ct level from the background 
region. Masked-out regions covered about 3.3% of the whole back- 
ground area. Finally we calculate the number counts from the re- 
main galaxies to estimate the background number counts in the 
cluster direction. 

Recent studies have provided ev idences that the LFs of 
GCs do vary with c lustercentric radius JSeiiersbergen et al]|2002l ; 
[Hansen et al] l2005l) . A suitable region should be chosen which 
is large enough to contain most member galaxies and do not 
include much projected contamination. We adopt a aperture of 
r=30'centered at the cluster centroid (about 1.7 Mpc at the rest 
frame of A671). The results of the background estimation and final 
luminosity function are showed in Figure 8. The apparent magni- 
tudes are converted to absolute magnitudes by the relation 



M = m- DM(z) - Ko.i{z) 



(4) 



Where DM{z) is the distance modulus as determined from the red- 
shift assuming a particular cosmology, i^o. i (z) is the /('-correction 
from a galaxy si z to z = 0.1. We esti mate the Ko.iiz) usin g 
the software KCORRECT(version 4-1-4. iBlanton et alj (12003')). 
As showed in the bottom-right panel of Figure 8, a single Schechter 
function can fit the data very well. The best fit parameters are 
0*=21.O, M*=-21.6, a =-1.12, which is in good agreement with 
Ide Filippis et al.l ( 1201 ih . No strong " upturn" at faint magnitudes is 
observed for A671. 



past decades an d it is well described by the Schechter function 
(ISchechteJ 19761) : 



<j3{L)dL = (ji* {L/L*)°'exp{-L/L*) d{L/L'), 



(2) 



where (j>*, L* , and a are the normalization parameter, the character- 
istic luminosity, and the slope parameter at faint end, respectively. 
In the domain of absolute magnitude, the Schechter function can be 
expressed as: 

,0.4(M*-M),^^ 



0(M)dAf = 0*lO''-*('^+^'(*^*-*^'ea;p[-lO° 



(3) 



where M* is the characteristic absolute magnitude. 

The most challenging issue in measuring the LFs of galaxy 
clusters (GCs) is that one needs to pick up cluster members out 
of background galaxies along the line of sight. Ideally, one needs 
the spectroscopic redshifts for all galaxies to exclude non-members 
in the cluster field. Unfortunately spectroscopic measurements are 
rather time-consuming. Our BATC photometry enables us to se- 



6 STAR FORMATION HISTORY AND ELEMENT 
ABUNDANCES 

As mentioned in the first section, many former studies sup- 
port that galaxy properties strongly correlate with local environ- 
ment. The remarkable morphology-segregation of A671 demon- 
strates that the morphologies of member galaxies strongly cor- 
relate with their local environment. In this section we would 
like to investigate the star formation histories (SFHs) of the con- 
firmed member galaxies by applying models of stellar population 
synthesis on their observed spectra. Two models will be taken 
in this section. On the one hand, we fit the SDSS spectra by 
the STARLIGHT code dCid Fernandes et alj|2005l : iMateus et all 
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Figure 9. Spectral synthesis of the brightest galaxies in A671. Top panel 
shows the obsei'ved spectrum (black line) and model spectrum (red Hne). 
Bottom panel shows the residual spectrum (black line) and the mask regions 
(pink lines). 



12004 ICid Femandes et alj|2007h to derive the physical parameters 
based on their SFHs. On the other hand, we compare the absorption 
line indices (the Lick/IDS indic es) of the member ETGs with the 
model predictions developed bv lShiavonl ( |2007|) (hereafter S07) to 
put constraints on the SFHs and chemical enrichment. 



6.1 Fitting Spectra with STARLIGHT 

We fit the SDSS spectra of member galaxies with the STARLIGHT 
codes, which aims at fitting an observed spectra with a liner combi- 
nation of theoretical simple stellar populations (SSPs). The model 
spectrum is given by 



N 



(g) G(u* 



(5) 



where Mx is the model spectrum, Mx„ is the synthesis flux at the 
normalization wavelength Aq, Xj is the so-called population vec- 
tor, bj^x is the jth SSP spectrum at A, and rx = 10"° '*<'*^"^^o) 
represents the reddening term. The G{v,,a,) is the line-of-sight 
stellar motions that modelled by a Gaussian distribution centered 
at velocity v, and with a dispersion of cr*. A'^ is the total number 
of SSP models. In our work, the SSP base is made up of A = 45 
SSPs, three metallicities (Z = 0.2^©, Zq, 2.5Zq) and 15 ages 
(from 1 Myr to 13 Gyr), which are taken from evolutionary mod- 
els in iBruzual & Chariot! (l2003l). The galactic extinction law of 
ICardelli. Clayton & MathisI ( Il989h with i?v = 3.1 is adopted. 

All SDSS observed spectra are shifted to the rest-frame, and 
then interpolated into a resolution of I A before fitting. Wave- 
length regions of emission lines are masked out. Figure 9 shows 
the spectral fitting for the brightest cluster galaxy (BCG) of A67I. 
As demonstrated by this figure, the combination of SSP spectra can 
fit the observed spectrum very well. 



The STARLIGHT presents the SSP fraction, intrinsic extinc- 
tion Ay, velocity d i spersi on a, and stellar mass M*. Following 
ICid Femandes e t aljjlool) . we derive the flux- and mass-weighted 
average ages, which are defined as 



< logt, >L= ^ajjlogfj; < logt, >^ 



JV 

E 



Ujlogtj, 



(6) 



where Xj is the flux-weighted population vector (i.e., the fraction 
of flux contributed by certain SSP), and uj is the mass-weighted 
population vector. The average metallicities <Zl> and <Zm> can 
be derived similarly. 

6.1.1 SFH via STARLIGHT Fitting 

As illustrated by ICid Femandes et aEI bOOS i). the individual out- 
put vector may be dramatically deviated from the simulated input 
value. However, the average values of stellar age and metallicity 
should be more reliable, whatever they are weighted by light or 
mass. The flux- weighted age is more sensitive to the young stel- 
lar component, so the mass-weighted age is more underlying and 
intrinsic. The situation is the same for average metallicity. In this 
section, we will take the average ages and metallicities weighted 
by stellar mass. Since stellar population in ETGs are dominated 
by the old components, the ETGs have their average ages within a 
relatively narrow range. 

Figure 10 presents the derived mass-weighted ages and metal- 
licities as functions of cluster-centric radius R and total stellar mass 
(M«) assembled in cluster galaxies. The upper two panels show 
the mass-weighted ages and metallicities as a function of R. The 
galaxies in the core region (R<400kpc) of A671 are denoted by 
filled symbols. The remarkable morphology-density segregation of 
A67I is well shown in these two panels. The LTGs (denote by blue 
symbols) locate in the outskirts (denoted by open symbols), and 
have younger stellar ages. On the other hand, the ETGs with older 
stellar ages are located in the core region. For the ETGs of A67I, 
no correlation is found between metallicities and R. 

The lower two panels of Figure 10 present how the mass- 
weighted ages and metallicities correlate with stellar mass. Both 
ages and metallicities are found to be correlated strongly with stel- 
lar mass, and such correlations are dependent upon morphology. 
In general, the more massive galaxies have older ages and richer 
metallicities. For the LTGs in A671, the linear correlations of age 
and metallicity with stellar mass appear tighter and steeper. Even 
for the ETGs with similar stellar mass, the ETGs in core region 
tend to have older ages than the outskirt ETGs. However, the ETG 
metallicities seem to not vary with cluster-centric radius. 

6.2 Lick Indices 

The age-metallicity degeneracy has haunted stellar population 
analysis for decades. Nevertheless, the promising approach to 
break i t remains the combi ned use of multiple absorption-line 
indices lKong & Cheng 2001). In this section, we will measure the 
absorption lines of ETGs in A67I and compare them with state- 
of-art SSP models in order to infer their ages, metallicities and a- 
enhancements. 

6.2.1 Lick Index Measurements 

The bandpass e s of the Lick indices are defined in Table I of 
IWorthev et al] h994l) . We measure the Lick indices with a mod- 
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Figure 10. Correlation of mass-weighted ages and metallicities witli cluster-centric radius R and stellar mass. Black circles represent the ETGs, while blue 
triangles represent the LTGs. Galaxies in the core region are denoted by filled symbols. 



ified ve rsion of the Lick_EW rout ine in EZ_Ages package devel- 
oped bv lGraves &Schiavonl(l2008l) Fl The Lick_EW routine reports 
the errors of each L ick index calculated in the way suggested by 
ICardiel et al] ( Il998h . The SDSS spectral resolution (69 km s"^) 
does not match the originally defined resolution of Lick indices. 
The Lick_EW routine smooths the SDSS spectra to the resolution 
of Lick/IDS system before measuring the indices. For the galaxies 
with high velocity dispersion, the smoothed absorption features in 
the SDSS spectra are at resolution even poorer than the Lick resolu- 
tion, the Lick_EW routine will apply cr-correction for these galax- 
ies. The output includes the measurements of cr-corrected indices 
and their errors. 

One of the most challenging issues in Lick index measure- 
ments is emission line contamination. The Balmer absorption fea- 
tures are contaminated by emission from ionized gas, either from 
star formation, AGN activity or interstellar shocks. Balmer line 
emission was estimated from equivalent of Ha, which are re- 
trieved from the MPA/JHU SDSS DR7 catalog. Some of previ- 
ous studies used E W(H/3)= 0.6 EW([0 in]A5007) for correction 
jTrag er et al. 2000). We do not use [OlIl]A5007 for correction be- 
cause Ne lan et al.. (,2005.) have found that the relation between H/3 
and [O in]A5007 for the ETGs varies in different mass ranges. They 
found that the correlation between H/? and Ha is much tighter. 
Emission EWs for higher order Balmer lines are obtained from 
Ha by assuming standard values from Balmer decrement (in the 

^ http ://w w w.ucolick. org/ ~ graves/EZ_Ages .html 



absence of reddening). In this way, EW(H/3)= 0.36 EW(Ha), 
EW(H7)= 0.19 EW(Ha), and EW(H5)= 0.13 EW(Ha). 



6.2.2 SSP Model and Stellar Population Parameters 

Our goal is to use the S07 model to derive the SSP-equivalent pa- 
rameters. Firstly we create grid to fit three parameters: ages, [Fe/H], 
and [Mg/Fe]. We use the solar-scaled isochrones and the Salpeter 
initial mass function suggested by the EZ_Ages documents. The 
[O/Fe] is set to be zero, and other a elements are tied to Mg. The 
ages range from 1.2 to 17.7 Gyr, and t he [Fe/H] ranges from -1.3 t o 
0.2. Other details can be found in S07. lGraves & Schiavoijjioosh . 

The three parameters are derived following two steps: i) 
for each galaxy, we firstly calculate age and [Fe/H] using H/3 
and <Fe>, at all [Mg/Fe]; ii) we compare the pair <Fe> 
(<Fe>=0.5(Fe5270-l-Fe5335)) and Mgb with models at the me- 
dian age obtained in the first step, and calculate a new [Fe/H] and 
[Mg/Fe]. We then update the ages by interpolating the ages found 
in the first step at the [Mg/Fe] derived in the second step, and it- 
erate the second step with the new age. Usually two iteration steps 
are needed before convergence. For those measurements beyond 
the model grids, we set the parameters to be the boundaries, i.e, the 
maximum or minimum of the models. 

Before deriving the stellar population parameters of our sam- 
ple with the SSP models, we compare our measurements to the pre- 
dictions of each model on the grids of age and metallicity. In Figure 
11, our measurements of H/3 and [MgFe]' are compared with the 
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[MgFe]' 

Figure 11. The index plot on the S07 model grid assuming [o/Fe]=0.2, 
shows the [MgFe]' vs. H/3. The filled and open circles represent the ETGs 
in the core region and those in outskirts, respectively. 



S07 models assuming [afPe]=0.2, where the index [MgFe]' is de- 
fined as follows: [MgFe]'=^Mgb(0.72Fe5270 + 0.28Fe5335). 
[MgFe]' is a good indicator of met allicity which is almo st indepen- 
dent upon a/Fe ration variations lf ThomasetaDl2003h (hereafter 
TMB). To convert between [Fe/H] and [Z/H], we adopt the rela- 
tion given by TMB: [Z/H]=[Fe/H]+0.94[a/Fe], and we assume 
[Mg/Fe]=[Q/Fe]. The median age is about 7 Gyr but with a large 
scatter. Most of the core-region galaxies have very small H/3 val- 
ues, and occupy the oldest end of age distribution. The [MgFe]' is 
mainly distributed along [Z/H]— 0, and does not have any environ- 
mental effects. 

We present the derived six SSP parameters as a function of 
velocity dispersion (a) in Figure 12. Galaxies which fall outside 
model boundary are not included in the figure. The parameter er- 
rors are given by the fitting procedure for each galaxy. We then 
compute an average error for each parameter by weighted the out- 
put error of individual galaxy with its signal-to-noise ratio. As pre- 
sented in Figure 11, majority of the galaxies falling out of the model 
boundary are located in the core, which makes hard to investigate 
environmental effect. Though only 46 ETGs are able to be fitted by 
the model, the relations between the SSP parameters and velocity 
dispersion (a) are still remarkable. 

In the top-left panel of Figure 12, the SSP ages show strong 
dependence upon velocity dispersion, and the low-cr galaxies span 
a wider age range, indicating that the low a galaxies have various 
possibilities of star formation histories compared with the high-a 
ones. Simi l ar resu lts are found in the bright galaxies in Coma by 
IPrice et al] ( 1201 ll) . It is noteworthy that the SSP ages derived by 
S07 model are not compatible with the average stellar ages given 
by the STARLIGHT fitting. Firstly, the SSP ages in STARLIGHT 
code range from 1 Myr to 13 Gyr, while the range of the SSP ages 
in S07 models is from 1 Gyr to 17.7 Gyr. The typical age of ETGs 
from STARLIGHT fitting (see Figure 10) is older than 7 Gyr. In 
Figure II, the ETG ages derived by the S07 model span a wider 
range. Additionally, the STARLIGHT gives the best-fit ages with- 
out errors. The S07 model determines the ages by comparing sev- 
eral combined-indices with the theoretical model, and the measure- 



ment errors can strongly affect the output ages. For the galaxies 
whose indices locate near the model boundary, measurement errors 
of indices will bring in greater uncertainties in the parameter fitting. 

The top-middle and top-right panels display the relation of 
[Fe/H] and [Mg/Fe] with a. The [Fe/H]— cr correlation is very 
tight, with a correlation coefficient of rg = 0.513, whereas the 
[Mg/Fe]— (T correlation is rather weak. The former studies have 
found [Mg/Fe] strongly correlated with a ('T homas et al] l2005l: 
IZhuetalj[2O10l) . We should keep in mind that the tightness and 
slope of linear correlation strongly depend on the sample size. Con- 
sidering our small sample size and the intrinsic scatter of this cor- 
relation, a relatively weaker [Mg/Fe]— cr relation is still reasonable. 

In the three bottom panels of Figure 12, we present [C/Fe]— cr, 
[N/Fe]— cr, and [Ca/Fe]— a relations. Only a weak correlation of 
[N/Fe]— cr is found, with a correlation coeffi cient of Vs = 0.404 . 
The results of our fitting are similar to those in lGraves et al. I( l2008h . 
and their results are based on a large sample of about 6000 red 
sequence galaxies from the SDSS. 



7 DISCUSSION 

We have investigated the dynamics of A67 1 based on the spectro- 
scopically confirmed members, and the result of K-test strongly 
suggests that A671 has significant substructures. Many authors 
have foun d that a l arge fr action of galaxy c lusters have sub- 
structures dPresslerl Il988l : iMohr et alj Il993l ; lYuan et al.l l2003l : 
I Yang et ^12004 ). indicating that massive clusters may assemble 
their masses and grow up by accreting small groups. The contour 
map of member g alaxies in A67I fits well with the X-ray inten- 
sity contour map. iRamella et ^ ((2003) suggested that there are 
two substructures in A67I. Location of one substructure is well 
associated with the potential substructure B, and the other is at the 
south part of cluster, which is not significant enough to be detected 
by the K-test. Their substructure-finding algorithm is based on the 
projected positions of galaxies, and does not utilize the redshift in- 
formation. 73% of the clusters in their sample were found to have 
substructures, and this fraction is higher than most studies. Their 
magnitude limit of galaxy samples is V ~21'".2, much fainter than 
our limiting magnitude /ibatc ~ 19™ .5, thus the projection effect 
cannot be ignored, and some substructures they found could be un- 
true. After 97 newly selected galaxies included, the substructure 
B is enhanced, while another substructure appears less prominent. 
The follow-up spectroscopy of these faint galaxies is need for re- 
vealing the details of dynamical substructures in A671. 

The remarkable morphology-density relation indicates that the 
cluster environment indeed have played an important role in evo- 
lution of cluster galaxies. A visual inspection of the C-M diagram 
shows that a large fraction of bright member galaxies (/ibatc < 
18™. 5) have evolved to be the "red sequence", and faint member 
galaxies with /ibatc ~ 19^.0 are found to have considerable ac- 
tivities of star formation. When the galaxies are accreted into a 
cluster, their star formation activities are expected to be suppressed 
by so me important proces ses, such as tidal stripping and "harass- 
ment" I Moore et_al.lll996h . ram pressure stripping of the gas disk 



( lAbadietalJl999l), and rem oval of gas reservoir surrounding each 
galaxy dBalogh et alj|2002h . 

We derived the average ages and metallicities for member 
galaxies by fitting their spectra. The most remarkable feature of the 
age distribution is that the E TGs in the cor e region have older ages 
than those in the outskirts. iThomas et all (|2005) have found that 
the ETGs in dense environment have average ages ~ 2 Gyr older 
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Figure 12. The .six SSP-equivalent parameters as the functions of velocity dispersion a for the 46 galaxies which can be fitted by the models. The symbols 
are defined in Figure 10. The black line represents the linear fit. 



than those in field environment. They derived ages by comparing 
the Lick indices with the prediction of TMB model. Our results 
confirm their conclusion though the average ages are derived by 
different methods. This can be interpreted in two ways. On the one 
hand, theoretical work shows that dark matter h alos in dense en - 
vironments were assembled earlier than average jGao et al ] l2005h . 
As a result, the galaxies in the core region of clusters formed earlier 
than those in the outer regions, and thus have older stellar ages. On 
the other hand, the older stellar age of the ETGs in the core region 
could be explained by lack of recent star formation compared with 
those in the outer region. Galaxy clusters have dense gas with high 
temperature, and the core-region galaxies lost their gas reservoirs 
by interacting with the dense intracluster medium (ICM) (usually 
by ram pressure stripping). Thus the core-region galaxies should be 
gas-poor, and have less possibilities of recent star-formation. Aside 
from age, the total metallicities of the ETGs show subtle dep en- 
dence on environment, which agrees well withl Zhu et alj jZOlOl) . A 
larger sample of cluster ETGs is needed to investigate the environ- 
mental effects on chemical evolution. 

The dependence of age upon stellar mass could be explain- 
able i n the downsizing scenario for galaxy formation jCowie et al.l 
Il996h . In this scenario, star formation lasts longer in less massive 
galaxies than in more massive galaxies. Thus massive systems will 
have older SSP ages on average. Evidences of downsizing effect in 
the local universe have been found by many recent studies on the 
stellar populati on of ETGs ( Nel an et al. 2005; Graves et al. 2008; 
IZhu et al.i20ld;|Price et alj201 ij). The age-(T slope in these studies 
span a range between 0.35 and 0.93. The different slope values are 
mainly due to the different sample property (e.g., galaxy type, sam- 
ple size, etc.) and methods of index correction (e.g., emission infill 
correction, velocity dispersion correction, etc.). Price et al. ( 2011) 
have tested the roubustness of their observed age-cr slope against 
these two factors when they studied the stellar population in Coma. 
In their test, a stricter emission-line cut and different methods of 
velocity dispersion correction were used. Only slight change on the 



age-cr slope was found in their tests. They concluded that their data 
robustly support the downsizing scenario. 

The ETGs also have very strong dependence of a-abundance 
on a. The a elements are mainly from the Type II supemovae, 
and the iron-peak elements come mainly from the Type la super- 
novae. The stronger alpha-enhancement in the more massive ellipti- 
cal galaxies may imply that their star formation timescale is shorter 
than less massive elliptical galaxies, before the delayed Type la su- 
pemovae enrich the star-forming regions with iron-peak elements. 
The observed a-abundance— cr relation fits well with the predictio n 
from the hierarchical models with feedback JPe Lucia et al ]|200^ . 



8 CONCLUSION 

This paper presents a photometric study of A671 with the Beijing- 
Arizona-Taiwan-Connecticut (BATC) multicolor system and the 
SDSS data. The main conclusions can be summarized as follow: 

(i) About 7000 sources are detected in a BATC field of 
58' X 58' centered at A671, and their SEDs in 15 intermediate 
bands are obtained. The 985 galaxies brighter /ibatc = 20^.0 
are selected by cross-identifying our BATC source catalog with 
the released catalog of SDSS galaxies. There are 205 galaxies with 
known spectroscopic redshifts in our viewing field, among which 
103 galaxies with 0.04 < z^p <0.06 are selected as spectroscop- 
ically confirmed members of A67 1 . The sample of bright member 
galaxies is composed of 63 ETGs and 40 LTGs. 

(ii) The dynamics of A671 is investigated based on the 103 
spectroscopically confirmed members. The result of K-test on dif- 
ferent scales strongly suggest that A67 1 have significant substruc- 
tures. Three potential substructures have been suggested with the 
method of localized deviation of velocity distribution. 

(iii) Photometric redshift technique is applied to the 985 galax- 
ies for further membership determination. Our photometric red- 
shifts (zph) of the bright members are basically consistent with the 
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spectroscopic redshifts (^sp)- Base on the statistics of photometric 
redshifts, the galaxies with 0.028 < Zph < 0.073 are selected as 
member candidates. After further selection by the color- magnitude 
relation, 97 galaxies down to ftsATC = 19™. 5 are picked up as 
faint members of A671. 

(iv) Based on the enlarged sample of member galaxies, spatial 
distribution and velocity structure of A671 are studied. Since the 
large Zph uncertainty of faint galaxies have smoothed the localized 
abnormity in velocity distribution, the K-test of the enlarged sam- 
ple does not confirm the three substructures mentioned above. The 
morphology-segregation becomes very remarkable after the faint 
members are taken into account. The luminosity function in the 
SDSS r-band shows a flat slope at faint end, a ~ —1.12. 

(v) Mass-weighted stellar ages and total metallicities of bright 
members are derived by fitting their spectra with the spectral syn- 
thesis code, STARLIGHT. The ETGs in the core region have older 
ages than those in the outskirts. The more massive ETGs are found 
to be older than the less massive ones. No environmental effect is 
found for the metallicities of the ETGs. Strong correlations of mean 
stellar age and metallicity with stellar mass are confirmed, and such 
correlations are found to be dependent upon morphology. The pos- 
sitive age-mass correlation supports the downsizing scenario. 

(vi) A set of Lick indices of the ETGs is measured in order 
to derived their SSP-equivalent stellar parameters (such as age, 
[Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe]) by utilizing S07 
model. The ETGs at cluster center tend to have smaller H/3 indices, 
indicating that central ETGs are likely to be older. The total metal- 
licity indicator [MgFe]' does not show any environmental effects. 
The relations between the six SSP-parameters and velocity disper- 
sion (fj) are also studied. The relations between the SSP-parameters 
and a in A671 are in good agreement with previous studies. 
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